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We study the thermodynamics of four-dimensional Kerr-Newman-AdS black 
holes both in the canonical and the grand-canonical ensemble. The stability 
Oh, conditions are investigated, and the complete phase diagrams are obtained, 

(— I \ which include the Hawking-Page phase transition in the grand-canonical en- 

semble. In the canonical case, one has a first order transition between small 



X 
J3 



and large black holes, which disappears for sufficiently large electric charge or 
angular momentum. This disappearance corresponds to a critical point in the 
phase diagram. Via the AdS/CFT conjecture, the obtained phase structure 
is also relevant for the corresponding conformal field theory living in a rotat- 
ing Einstein universe, in the presence of a global background U{\) current. 
An interesting limit arises when the black holes preserve some supersymme- 
try. These BPS black holes correspond to highly degenerate zero temperature 
states in the dual CFT, which lives in an Einstein universe rotating with the 
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I. INTRODUCTION 



The conjectured equivalence of string theory on anti-de Sitter (AdS) spaces (times some 
compact manifold) and certain superconformal gauge theories living on the boundary of AdS 
has lead to an increasing interest in asymptotically anti-de Sitter black holes. This 
interest is mainly based on the fact that the classical supergravity black hole solution can 
furnish important informations on the dual gauge theory in the large N limit, denoting 
the rank of the gauge group. The standard example of this is the well-known Schwarzschild- 
AdS black hole, whose thermodynamics was studied 17 years ago by Hawking and Page [^, 
who discovered a phase transition from thermal AdS space to a black hole phase, as the 
temperature increases. This means that at a certain temperature thermal radiation in AdS 
space becomes unstable, and eventually collapses to form a black hole. The Hawking-Page 
phase transition was then reconsidered by Witten in the spirit of the AdS/CFT corre- 
spondence. He observed that it can be interpreted as a transition from a low-temperature 
confining to a high temperature deconfining phase in the dual field theory. 
Up to now, the study of black holes in the context of the AdS/CFT correspondence has been 
extended in various directions, e. g. to black holes with non-spherical event horizon topolo- 
gies or to electrically charged Reissner- Nordstrom- AdS (RNAdS) solutions |lTO| , pA| , 
following previous work on this subject performed in Thereby, some interesting ob- 



servations have been made, like the striking resemblance of the RNAdS phase structure to 
that of the van der Waals-Maxwell liquid-gas system [|TU],|rT|[|, or the possible appearance of 
so-called "precursor" states [|1^] in the CFT dual of hyperbolic black holes 0. Furthermore, 
the thermodynamics of R-charged black holes in four, five and seven dimensions has been 
studied in detail in [^]. Another interesting extension is the inclusion of rotation, i. e. the 
consideration of Kerr- AdS black holes [lT8| , |l9| pl. In this case, the dual CFT lives in a rotating 



Einstein universe |T^. The authors of |T8| studied the limit of critical angular velocity, at 
which this Einstein universe rotates with the speed of light, from both the CFT and the 
bulk side. Recently, the AdS/CFT correspondence has been probed in |21]] by comparing 
in more detail the thermodynamics of rotating black holes in five-dimensional anti-de Sitter 
space with that of A/" = 4 super Yang-Mills theory on a rotating Einstein universe in four 
dimensions. 

In the present paper, we shall consider the charged rotating case, i. e. Kerr-Newman-AdS 
(KNAdS) solutions. (Note that the appearance of second-order phase transitions for Kerr- 
Newman black holes in de Sitter space, i. e. for positive cosmological constant, was first 
shown by Davies [p^] ). 

The remainder of this article is organized as follows: 

In section |T| we review the four-dimensional KNAdS black hole, give the conserved quanti- 
ties like mass and angular momentum, and recall the conditions under which some amount 



similar phase structure has recently been found by considering stringy corrections to AdS black 
holes H. 



For a discussion of Kerr-NUT-AdS and Kerr-bolt-AdS spacetimes cf. [20|. 
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of supersymmetry is preserved. In |T| we calculate the Euclidean action, using the method 
of counterterms introduced in |^-26]. In section IV the black hole thermodynamics is con- 
sidered in both the canonical and grand canonical ensemble. The resulting phase structure 
is elaborated, and the thermodynamical stability of the solutions is studied. Finally, in |^ 
our results are summarized and discussed in the context of the conformal field theory in the 
rotating Einstein universe. 



II. KERR-NEWMAN-ADS BLACK HOLES 



Here we consider the usual charged rotating black hole in AdS space. Its horizon is 
homeomorphic to a sphere, and its metric, which is axisymmetric, reads in Boyer-Lindquist- 
type coordinates P7| , p5| 



dt - 



asm 9 



+ ^dr' + ^d9' + 



Ao sin^ e 



a dt 



2 I 2 H 2 



where 



2 2,2 2 n 

p = r + a cos 



(2) 



A, = (r^ + a^) 1 + 



P 



— 2mr + z 



(3) 



Here a denotes the rotational parameter and z is defined by z'^ = + 1e Qm being 
electric and magnetic charge parameters respectively. 

The metric ([^) solves the Einstein-Maxwell field equations with an electromagnetic vector 
potential and an associated field strength tensor respectively given by 



and 



A 



QeT 1m COS 9 3 

e ^= e , 



PV A0 sin 9 



F = J [ge (r ^ - cos^ 9) + 2g^ra cos 9] e° A e^ 



(4) 



P 



where we have defined the vierbein field 



H — - [gm(r^ - cos^ 9) - 2qera cos 9] e^ A e^ 



^A^ , , a sin^ 9 , 
-(dt ^d0). 



P 



/a: 



dr, 



(5) 



(6) 



d^, 



' Ag sin 9 
1 

P 



a dt 



+ 



(7) 
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Let us define the critical mass parameter niextr, 




i + T^j +7^(«' + ^')-7j-i| • (8) 



A study of tlie positive zeroes of tlie function A,, shows that the line element (|I|) describes a 
naked singularity for m < rUextr^ and a black hole with an outer event horizon and an inner 
Cauchy horizon for m > m^xtr- Finally, for m = mextri the lapse function has a double root 
and represents an extremal black hole. One should also note that the solution is valid 
only for a? < it becomes singular in the limit a? = l"^ . In this critical limit, which has 
been studied extensively in Ref. [0] (cf. also |^), the three dimensional Einstein universe 
at infinity rotates with the speed of light. Here we assume m to be larger than rrtextr and 
a < Z, so the metric (p represents an AdS black hole, with an event horizon at r = r+, 
where r+ is the largest solution of = 0. For the horizon area one gets 

A = l^fiLti!!). (9) 

Analytical continuation of the Lorentzian metric by t ^ ir and a ^ ia yields the Euclidean 
section, whose regularity at r = r+ requires that we must identify r ~ r + /9 and ~ 
+ ifi^Hi where the inverse Hawking temperature (3 is given by 

, "-^^f^ .... . (10) 



l + jT + Ajr 



and 



r| + 

represents the angular velocity of the event horizon. ^2^;^ can also be obtained by using the 
fact that the surface r = r+ is a bolt of the co-rotating Killing vector dr + iVLhO^ in the 
Euclidean section. One can also write the metric ([^) in the canonical form 

ds^ = -NHe + ^dr^ + Cde^ + ^L^^{d<P - udt)\ (12) 
A^ Ae p^^^ 

where 

S2 = (r2 + a2)2Ae-a^A,sin2^, (13) 



^Note that the solution with m = z = describes AdS space seen by a rotating observer. 



5 



and the lapse function and the angular velocity uj are defined by 



= (14) 



UJ 



[Ae{r' + a^) - Ar] (15) 



respectively. Note that uj = Qh on the horizon, whereas u = —a/t^ = Q^o for r oo. The 
fact that the angular velocity does not vanish at infinity is a salient feature of rotating black 
holes in AdS space, different from the asymptotically flat case, where floo = 0. The angular 
velocity fl entering the thermodynamics is then given by the difference 

n = nH-n^, (i6) 



analogously to the case of an electric potential. For asymptotically flat spacetimes, (|I6| 
reduces to = Qh, but for the KNAdS black hole we get 



This is exactly the angular velocity of the rotating Einstein universe at infinity ITsf. To see 



this, one first applies an implicit coordinate transformation, which takes the standard AdS 
metric to the m = z = KNAdS form [Q. ([T7|) can then be read off from the transformation 



rule of the coordinate (p. 

The fact that the angular velocity relevant to Kerr-Newman-AdS black hole thermodynamics 
turns out to be the one of the rotating Einstein universe at the AdS boundary, agrees nicely 
with the AdS/CFT correspondence: If the KNAdS black hole in the bulk is described by a 
conformal field theory living on the boundary, then the relevant angular velocity entering 
the thermodynamics should be that of the rotating Einstein universe at infinity. 
Moreover, when Q < l/l, a. timelike killing vector can be globally defined outside the event 
horizon , in contrast to the asymptotically flat case. Then the rotating black hole can be 
in equilibrium with rotating thermal radiation all the way out to infinity, and a consistent 
thermodynamics can be defined. When this condition is fullfilled, the rotating Einstein 
universe spins slower than light, and an associated state can be defined on the boundary 
CFT. This is closely related with the fact that for Q < l/l there are no superradiant modes 

H. 



The mass M and the angular momentum J can be defined by means of Komar integrals, 
using the Killing vectors dt/r. and d^^. Taking AdS space as reference background, one gets 

M = — , J=i^- (18) 

The charges are obtained by computing the flux of the electromagnetic field tensor at infinity, 
yielding 



^We have normalized the Killing vectors so that the corresponding conserved quantities generate 



the algebra so(3, 2), like it was done in Ref. [29|. 
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The electric potential measured at infinity with respect to the horizon, is defined by 

$ = KxX^^ - Alx\=r, = -1^, (20) 

where x = 9t + ^ndfp is the null generator of the horizon and denotes the electric part 
of the vector potential. 

Finally, let us recall the conditions under which some amount of supersymmetry is preserved 
by the Kerr-Newman-AdS black holes. These objects were considered in p9| , p0| in the 
context of gauged N = 2, D = 4 supergravity. There it was found that for 

m^ = la(l + jy, q^^=la(l + jy, = (21) 

the KNAdS black holes preserve half of the supersymmetries. In the BPS limit we have 



via ^ via ^ ^ ylaa 

M=— Qe = ^, Qm = 0, J=--—,, (22) 

(1 - f ) 1 - T (1 - f ) 



so the Bogomol'nyi bound is |^ 



M = Qe + y. (23) 

Furthermore, in the supersymmetric case, one has the relations 

rl = al, Q = J. (24) 

The latter relation means that in the BPS limit, the Einstein universe at infinity, in which 
the dual conformal field theory lives, rotates effectively with the speed of light. An analogous 
situation was found in Ref. [ITS] for the supersymmetric rotating BTZ black hole. 



III. EUCLIDEAN ACTION 

We turn now to the calculation of the Euclidean action, which will then yield the ther- 
modynamic potentials relevant to the various ensembles In the action computation 
one usually encounters infrared divergences, which are regularized by subtracting a suitably 
chosen background. Such a procedure, however, in general is not unique; in some cases the 



choice of reference background is ambiguous, e. g. for hyperbolic AdS black holes p2| , pl 
Recently, in order to regularize divergent integrals like those appearing in the computation 
of the Euclidean action, a different procedure has been proposed PB|-PB|. This technique 
was inspired by the AdS/CFT correspondence, and consists in adding suitable counterterms 
let to the action. These counterterms are built up with curvature invariants of a boundary 
dAi (which is sent to infinity after the integration) , and thus obviously do not alter the bulk 
equations of motion. This kind of procedure, which will also be employed in the present 
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paper, has the advantage to be free of ambiguities, which, on the contrary, are present in 
the traditional approach in some particular cases, like mentioned above. 

To start with, we write the Euclidean Einstein-Maxwell action in the form 



^bulk ~\~ I surf ~\~ let 



^— [ d'x^[R-2A-F'] 

--^ / d^xVhK + let , 
^'^(^ JdM 



(25) 



A = —3/P being the 4-dimensional cosmological constant, with / the radius of AdS space. 
The action (|25| ) differs from the familiar one by the presence of the last term, which contains 
all surface counterterms needed to assure convergence of the integrals. In four dimensions 
it reads \M\ 



I'Ct 



dM 



I 



I' 



- + -n 

I 2 2 



ab 



(26) 



hab and TZab denoting the metric and Ricci curvature of the (arbitrary) boundary dM.. 
The variation of the bulk integral with respect to the fields gives 



61 



Gab + ^9a 



"^Fa^Fhc — -F gab 



^9 



ab 



, ^ / d'x^ {VaF^') SA, + -^ [ d'x Vh UaF'^'SA, 
47rG 47rG Jgj^ 



(27) 



where is the outward pointing unit normal to dAi. Asking for a stationary point of the 
action, the bulk integrals in Eq. (^) yield the Einstein-Maxwell field equations 



Gab + Agab = Tab , aF"" = 0, 

where we have defined the electromagnetic stress tensor 

Tab = 2F'^aFcb — -F^gab , 



(28) 



(29) 



while the surface integral has to vanish to have a differentiable action functional and a well 
defined action principle. This imposes the boundary condition 6Aa = on dAi. This action 
is thus appropriate to study the ensemble with fixed electric potential $e and fixed magnetic 
charge Qm- 

To study the canonical ensemble with fixed magnetic and electric charge, we have to add 
another boundary term to impose fixed Qe as a boundary condition at infinity , restoring 
the electromagnetic duality. The appropriate action in this case is 



/ = / - / d^x Vh HaF'^'Ab 

and its variation reads 

SI = (bulk terms) 



1 



AttG 



d^x Vh A^S (n^K 



ab 



(30) 



(31) 



dM 
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The bulk contribution yields again the Einstein-Maxwell equations, while the vanishing of 
the surface term, needed to have a differentiable action functional, requires 6{n°'Fab) = at 
infinity boundary condition. 



Now the evaluation of the Euclidean action is a straightforward computation. First of 
all, on shell the bulk contribution simplifies to 



Ibulk 



16nG 



(fx ^ 



(32) 



Then one chooses the boundary dJ^ = x S'^, S'^ being a 2-sphere with a large radius, 
which has to be sent to infinity after the integration. Since the metric is stationary, the 
time integration gives rise to a simple multiplicative factor (3. The integration on the other 
variables requires a little bit of work, but anyway it can be performed in a closed manner. 
The final results assume the form 



/3 



-ri + + 



+ 



(33) 



/3 



4GI 



2^ 



-ri + + + 2 



(34) 



valid for fixed potential and fixed charge respectively. In the following, we will consider only 
the case of vanishing magnetic charge, qm = 0, so we set qe = q and Qe = Q- 

The behaviour of the obtained actions as functions of r_|_ determines the thermodynamical 
properties of the black holes. It will be discussed in detail in section Note that both 
actions are singular in the extremal case (m = rriextr) and also for S = 0. The singularity 
for extremal black holes comes from the divergence of the inverse temperature jS. Since one 
divides the action by (3 in order to obtain the thermodynamical potentials, the latter are 
well-defined also for T = 0. 

In the last section we have already indicated how to calculate physical conserved quantities 
by evaluating Komar integrals. This traditional technique involves the subtraction of a 
suitably chosen reference background (in our case AdS space), in order to regularize infrared- 
divergent integrals. In the approach we are considering, instead, one does not use any 
reference background, and so it seems quite natural to compute the conserved charges by 
the direct use of the finite action (EBI), applying the method developed by Brown and York 



Let us briefly describe how this formalism works. We indicate by u"" the unit normal 
vector of a spacelike hyper-surface ^St at constant t and by dAi the spatial boundary of 
the spacetime manifold Ai, with induced metric hat- Moreover, S = 5*^ is the spacelike 
intersection ^St fl dA4 embedded in dj^ with induced metric aab- 

One starts by deriving the local surface energy-momentum stress tensor 
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which characterizes the entire system, including contributions from gravitation and (if 
present) matter. It is related to the stress energy tensor T"^ = -^-^-^ in a rather compli- 
cated way. 

Then, for any Killing vector field associated with an isometry of the boundary three- 
metric, one defines the conserved charge Q 



(36) 



Now, a straightforward application of (35) with the action (25) gives 

1 



SttG 



(37) 



The first term on the right-hand side of Eq. ( P7D results from the boundary term in the 
action, while all other terms are due to the presence of the counterterms we added in order 
to have finite quantities, when the boundary is sent to infinity. Employing the AdS/CFT 
correspondence, the result (|37|) can also be interpreted as the expectation value of the 
stress tensor in the boundary conformal field theory. In the case of the Kerr-Newman-AdS 
geometry, we choose dM. to be a three-surface of fixed r, and obtain for 

SnGru = ^ + Oiiy) , 
rl 



2ma 



rl 
m 



StiGtss = sin^ + sin^ 9 - a^] 

SnGree = + ^(l/r^) , 
rAo 



C(l/r2) 



(38) 



all other components vanishing. As expected, the presence of counterterms cures also the 
divergences of the charges. In fact we get 



m 



M = Q9,/s = - 



(39) 



T-n - """^ 



(40) 



in full agreement with the results (18) obtained with Komar integrals. 



^Note that in the conventions of [p3] an additional minus sign appears in the conserved charges 



(|3q); this comes from the fact that the definition of the extrinsic curvatures in |34] differs from 
ours by a factor of —1. 
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IV. THERMODYNAMICS 



A. Generalized Smarr Formula 



Using the expressions ([T^) for mass and angular momentum, for the horizon area, and 
the fact that = for r = r^, one obtains by simple algebraic manipulations a generalized 
Smarr formula for Kerr-Newman-AdS black holes, which reads 



A 



J2 



levr A 



+ -{AJ' + Q^) + ^ + — + 



A 



A_ 



A' 



327r2/2 



(41) 



In the limit / —>■ oo this formula reduces to the usual Smarr formula for asymptotically flat 
Kerr-Newman solutions [^,^, so one can consider the last two terms on the right-hand 
side as AdS corrections. 

Relation (^Tj) holds for classical black holes, and contains as usual all the information about 
the thermodynamic state of the black hole, taking the black hole entropy to be one quarter 
of the horizon areaPI, 



(42) 

and 



In the next subsection, we shall check Eq. (^2|) using the Euclidean actions (^), (| 
standard thermodynamical relations. 

We may then regard the parameters S, J and Q as a complete set of energetic extensive 
parameters for the black hole thermodynamical fundamental relation M = M{S, J, Q), 



S 



TT 



(43) 



One can then define the quantities conjugate to S, J and Q. These are the temperature 



T 



dM 



dS 



JQ 



8nM 



the angular velocity 



dM 



dJ 



SQ 



7lJ 



_s_ 



(44) 



(45) 



and the electric potential 



dM 



dQ 



SJ 



2MS 



Q q2 



(46) 



These expressions reduce to the corresponding asymptotically flat Kerr-Newman expressions 



in the / oo limit |36]. The obtained quantities satisfy the first law of thermodynamics 



^In the following we set G = 1. 
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dM = TdS + ndJ + $dg. 



(47) 



Furthermore, by eliminating M from (^)-(^) using (^), it is possible to obtain three 
equations of state for the KNAdS black holes. 

One easily verifies that the relations (^i]) , (^51) and (|46|) for temperature, angular velocity and 
electric potential respectively, coincide with equations (0), (p^ and ( pO| ) found in section 

Another quantity of interest is the thermal capacity Cj^q at constant angular momentum 
and charge, which, as we shall see in section |IV Q is relevant in the stability analysis of the 
canonical ensemble. It reads 



C, 



T 



dS 
df 



A-R-^MTS^ 



JQ 



A-kMTS 
1 - AtiT (2M + TS) + ^ (g2 + 



6S2 

1^ 



(48) 



In the above discussion, we have treated the cosmological constant as a fixed parameter; 
however, as shown by Henneaux and Teitelboim [R^, it is possible to induce it from a three 
form gauge potential coupled to the gravitational field. The interest of this mechanism 
resides in that it follows from Kaluza-Klein reduction of supergravity theories [^, and in 
particular it is relevant for the compactification of M-theory on S*^. We can thus promote 
the cosmological constant to a thermodynamic state variable. Its conjugate variable reads 



e 



dM 



dA 



SJQ 



2M 



3 



S 



S 



TT 



Ag3 
187r3 



(49) 



Now we can consider ensembles where A is allowed to fluctuate, and the first law reads 

dM = TdS + fidJ + $dQ + OdA. (50) 

Finally, if we regard M as a function of 5, J, and A, it is a homogeneous function of 
degree 1/2. Applying Euler's theorem we obtain 



hi 

2 



TS + nj+-^Q 
TS + nj + ]^<^Q 



BA 

s 

2^ 



$5 nj 
^^iT^ 



(51) 



Again, sending / to infinity we recover the usual Smarr law, and the last term is an AdS 
correction. 



B. Thermodynamic Potentials from the EucHdean Action 

We now turn to the definition of the thermodynamic potentials from the Euclidean ac- 
tions (|3^) and (p3D, relevant for the canonical and the grand-canonical ensemble respectively. 
Let us first treat the latter case. The Gibbs potential G(T, ^2, $) is defined by 



12 



GiT,n,^) = -, 



(52) 



with / given by (|33|). Using the expressions (|T^) for the temperature, ( p!7D for the angular 
velocity, and (|20|) for the electric potential, we get after some algebra the corresponding 
extensive quantities 



S = - 



dG 
df 



J = - 



n<s> 



dG 



Q = - 



dG 
9$ 



(53) 



TQ 



which turn out to coincide precisely with the expressions (^2]), (|T^) and (pJ]). One further 
readily verifies that 

G(r, n, ^) = M -nj -TS - <i>g, (54) 

which means that G is indeed the Legendre transform of the energy M{S, J,Q) (cf. (|43|)) 
with respect to S, J and Q. 

As for the canonical ensemble, the Helmholtz free energy is defined by 



F{T,J,Q) = - + nj, 



(55) 



with I given by (p^ . The term fij comes from the needed Legendre transformation to 
hold the angular momentum fixed. Using the relations ([T0| ) for the temperature, (|T8|) for 
the angular momentum, and (|19D for the electric charge, one verifies that the conjugate 
quantities 



dF 

^~ dT 



JQ 



dJ 



TQ 



$ = — 
dQ 



TJ 



agree with expressions (^2|), (|T^) and (|20|) . Furthermore one easily shows that 

F(T, J,Q) = M- TS, 
so F is in fact the Legendre transform of M{S, J, Q) with respect to S. 



(56) 



(57) 



C. Stability Analysis 

The stability of a thermodynamical system with respect to small variations of the 
thermodynamic coordinates, is usually studied by analyzing the behaviour of the entropy 
S'(M, J,Q). The stability requires that the entropy hypersurface lies everywhere below its 
family of tangent hyperplanes, as to say that the entropy must be a concave function of 
the entropic extensive parameters. This will place some restrictions on physical observables. 
For example, the thermal capacity must be positive in any stable system. The stability can 
also be studied by using other thermodynamic potentials, as for example the energy or its 
Legendre transforms, which have to be convex functions of their extensive variables, and 
concave functions of their intensive variables. The use of one thermodynamic potential with 
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respect to another is of course a matter of convenience, depending on the ensemble one is 
deahng with. Here we are mainly interested in finding the zones where the system is locally 
stable. These are bounded by "critical" hypersurfaces, on which 

det (^^irr) = 0, (58) 



dXidXj 
where Xi = M,J,Q. 

In the canonical ensemble, charges and angular momentum are fixed parameters and for this 
reason, the positivity of the thermal capacity Cj^q is sufficient to assure stability Jl^- This 



means that the points (hypersurfaces) where Cj^q vanishes or diverges represent the critical 
hypersurfaces we are looking for. The equation of such surfaces can easily be obtained by 
deriving the expression for T, Eq. (|^) (where M is given by (^)), with respect to S and 
requiring it to vanish (or to diverge). One can see that the equation dT/dS = 0, oo has only 
one physically acceptable solution, which reads [] 

2 _ _3S\7T + S)\7f + 2S) + 7rV(3vr + 2S) + 2t^^Q^S{2t^ + 35)) 
~ 47r4(37r + AS) 

^ ^ -K^Q^ + Q-K^Q^SiTT + 5)' - 27r^Q6(27r + 3S) 



2tt^{3tt + 45) 

+S\7T + 5)^(37r3 + IOtt^S + 157fS^ + 9S^) 

1 

+ 7r^g^(47r^ + 977^5 + Gn^S^ + Ott^^ + 6S^) ' . (59) 

In the grand-canonical ensemble the critical surfaces can be determined by requiring that the 
determinant of the whole Hessian of the Gibbs potential G{T, Q, $) be vanishing (diverging) 



17| . This yields an equation which can easily be resolved for the charge parameter q. The 



result is 

2 _ -a^ -2a^rl-2a*rl-rl-a*rl + 2rl + 2a'^rl + 3rl 
^ " -3a^ + rl + a^rl + ri ' 

These stability conditions will be analyzed in the following subsections. 



D. Canonical Ensemble 

In the canonical ensemble, we study the black holes holding the temperature T, the 
angular momentum J and the charge Q fixed. The associated thermodynamic potential is 
the Helmholtz free energy F{T,J,Q) (^5]). We shall analyze the thermodynamics in this 
ensemble in the (T, J)-plane, keeping the charge fixed. It immediately follows that in this 
plane the extremal black hole solutions are represented by the ordinate axis T = 0, and the 
half-plane T > corresponds to the nonextremal solutions. 



^In the following we set / = 1. 
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We shall begin by considering the local stability of the black holes. As explained in Section 



IV Q , local stability in the canonical ensemble is equivalent to the positivity of the thermal 
capacity (PS]). The behaviour of the thermal capacity can be more easily understood from 
the state equation T = T{S, J, Q) at fixed J and Q, Eq. In figure |I] we have reported 




FIG. 1. State equation T = T{S). The plot shows the T{S) curves at fixed Q = for J = 
(the upper curve), 0.005, 0.01, 0.0236, 0.05 and 0.1 (the lower curve). For « 0.0236 the local 
extrema merge in a point of inflection, and for higher J the curve is monotonically increasing. For 
Q / 0, see the text. 

these curves for Q = and different values of J. When J = 0, the usual Schwarzschild- 
AdS behaviour is reproduced; the curve first decrases towards a minimum, corresponding 
to the branch of small unstable black holes, then increases along the branch of large stable 
black holes. As soon as J 7^ 0, a branch of stable small black holes appears, separated 
from the branch of large black holes by intermediate unstable black holes. This results in 
the appearance of two phases in the canonical ensemble, a small black hole phase and a 
large black hole phase, in analogy to what was found in [10| for charged AdS black holes 
and in p[ for stringy-corrected black holes. While J grows, the local maximum of ^(5*) 
decreases, and eventually degenerates for Jc ~ 0.0236 into a point of inflection, signalling a 
second order phase transition for the Kerr-AdS black hole. For J > the curve T{S) is 
monotonically increasing, and we are left with a unique phase of black holes. The behaviour 
near this critical point is completely analogous to a liquid/vapour system described by the 
Van der Waals equation. 

In the Kerr-Newman-AdS case, one can repeat this analysis along the same lines; for < 
Q < Qc, where Qc ~ 0.166, we obtain a behaviour analogous to that of the Kerr-AdS black 
hole, with the difference that the stable small black hole branch appears already at zero 
angular momentum [|1^]. It is then possible to obtain the critical value JdQ), where the 
point of inflection is located, as a function of the charge. For Q = Qc the small black hole 
phase disappears, and for larger values of the charge the function T{S) is monotonically 
increasing for any angular momentum, and hence a unique stable black hole phase exists; 
the curve of second order critical points ends in J = 0, Q = Qc- The resulting phase diagram 
is shown in figure Q the solid line corresponds to the second order critical points. Inside 
the curve we have a region where two phases of black holes are allowed, with a first order 
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FIG. 2. Phase Diagram in the (Q, J)-plane and stability diagram in the (fi, J)-plane. In the left 
plot the curve corresponds to a line of second order critical points. Outside the region bounded 
by it, there is a unique black hole phase. Inside the region, one has a small black hole phase at 
low temperatures and a large black hole phase at high temperatures, separated by a first order 
phase transition. In the right plot we have drawn the lines where the thermal capacity diverges, 
for Q = 0, 0.05, 0.091, 0.12, 0.15 and 0.16 respectively from the outer to the inner curve. Outside 
these regions the thermal capacity is positive. The vertical line represents the limit where = 1/2; 
for Q > 0.091 the entire instability region has fl < l/l. 

transition between a small black hole phase and a large black hole phase as the temperature 
is increased, the two phases being separated by a coexistence curve. 

In the second graph of figure ^ we have visualized the instability region in the {fl, J)-plane. 
In the uncharged case we have stable black holes for large J and thermodynamic instability 
below the curve. As Q grows, the instability region shrinks and disappears for the critical 
charge Qc- The points inside the instability region correspond to mixed phases of small/large 
black holes, obtained by means of a Maxwell-like construction. 

The phase diagrams in the (T, J)-plane for the canonical ensemble are shown in figure 
^. The left diagram corresponds to the uncharged Kerr-AdS black hole, while the one on 
the right side is plotted for Q = 0.05 and shows the qualitative behaviour of the Kerr- 
Newman-AdS black holes with electric charge ranging in the interval < Q < Qc- In 
these diagrams, we have plotted the curves on which the specific heat diverges. Outside the 
region bounded by them, every point (T, J, Q) of the diagram corresponds to a single black 
hole solution, which is locally stable (its thermal capacity is positive). In contrast, inside 
the regions bounded by these curves, which are indicated by C < in the figure, every 
point (T, J, Q) corresponds in fact to three different black hole solutions. To elucidate this 
point let us go back to the curves T{S) of figure |I|. The divergence of the thermal capacity 
obviously corresponds to the local extrema of T{S). For every temperature between that of 
the local minimum and that of the local maximum of T{S), we have three black holes: a 
small metastable black hole, an unstable intermediate one, and a large stable one, which are 
represented by a single point in the (T, J, Q)-diagram. If we start from a small stable black 
hole and increase its temperature while holding J and Q fixed, as we reach the temperature 
of the local minimum of T{S), it undergoes a first order phase transition with a jump in the 
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FIG. 3. Phase diagrams in tlie {T,J)— plane. The phase diagram of the Q = Kerr-AdS 
black hole is plotted on the left, while the diagram for Kerr-Newman-AdS black holes with charge 
< Q < Qc ^ 0.166 is shown on the right. For higher charge the diagram is trivial, as one is 
left with a unique stable black hole phase. The left branch of the solid curve is the first order 
small/large phase transition curve. It ends in a cusp, where a second order phase transition occurs. 



entropy, and reaches the black hole represented by the local minimum of T{S). Then, the 
state will follow the branch of large stable black holes. Hence the low temperature branch 
of the CjQ = oo curve corresponds to a first order phase transition line between a small 
black hole phase and a large black hole phase. This line ends in a critical point located in 
the cusp, where the transition degenerates to a second order phase transition. Then, inside 
the region bounded by the Cjq = oo curve, the thermodynamically stable state is simply a 
locally stable large black hole phase. Finally, for Kerr-Newman-AdS black holes with charge 
Q > Qc the instability region and the small/large transition disappear, leaving a unique 
phase of stable black holes. 

The local stability is however not sufficient to ensure global stability; the full problem is 
very hard to tackle, but we could at least try to compare the free energy of the black 
hole with that of some reference background like thermal AdS space, which can lead to a 
Hawking-Page phase transition [Q. In doing so, however, one encounters the problem that 
pure AdS space is not a solution of the Einstein-Maxwell equations with fixed electric charge 
Q as boundary condition. Furthermore, AdS space has no angular momentum, as J = 
for m = (cf. ([ISD ). Therefore, it is not a suitable reference background in the canonical 
ensemble, which has fixed Q and J. One should compare the free energy of the black hole 
with that of other possible solutions, such as AdS space filled with a gas of charged particles 
carrying also angular momentum J, or a Kerr-Newman-AdS black hole with part of the 
total charge and angular momentum carried by a gas of such particles gravitating outside 
its horizon |]lT|. Only ioi Q = J = 0, i. e. the Schwarzschild-AdS black hole, pure AdS space 
contributes to the path integral and should be taken into account, as it was done in [Q. In 
the charged rotating case, however, the question of global stability and Hawking-Page phase 
transitions remains unsettled. 

Finally, we have to examine the stationarity condition Q < l/l, which is needed to have 
thermodynamic equilibrium, and a well-defined holographic thermal state in the boundary 
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CFT. The curve Q = l/l can be defined parametrically in the (T, J)-plane by the equations 
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with 



Se 



-(yr+ig^-i 



-oo 



(61) 



(62) 



The lower bound corresponds to supersymmetric black holes for Q ^ 0, and for S* — >■ oo 
the curve is always asymptotic to T = l/(27r/). The behaviour of the curve depends on 




FIG. 4. Q = l/l curves. These curves have been plotted for the values Q = 0, Q = 0.2 and 
Q = 0.5 of the charge, from the right to the left respectively. The vertical line is the asymptote 
T = l/(27rZ). For charges < Q < l/\/8 the qualitative behaviour is similar to the intermediate 
curve, while for higher charges the corresponding curves behave like the one with Q = 0.5. In 
region / we have Q > l/l and in region // we have il. < l/l. The intersection of these curves with 
the T = axis corresponds to supersymmetric black holes. 

the electric charge, and has been summarized in figure ^ For Q = 0, the curve is strictly 
decreasing from the asymptote and ends at J = for infinite temperature. If < Q < l/V^, 
the infinite temperature branch is deformed, curves back towards the J axis and intersects 
it in the supersymmetric point. Finally, for Q > 1/v^, the curve remains at T < l/27r/, 
and is strictly increasing from its intersection with the ordinate axis, again describing a 
supersymmetric black hole, to the asymptote. In any case the Q = l/l curve splits the 
(T, J)-plane into two regions, named / and // in figure ^. Black holes in region /, at the 
left of the Q = l/l curve, spin at angular velocity fl greater than l/l and show superradiant 
instability. On the other hand, black holes in region //, at the right of the curve, have 
Q < l/l, and have a well-defined thermodynamics. However, this is not the whole story. 
Inside the regions of local instability shown in figure |^, every point (T, J, Q) corresponds in 
fact to three black hole solutions: an unstable, a metastable and a locally stable black hole. 
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It is possible to verify that the locally stable one spins at Q < l/l, and hence the whole 
region of Q < l/l black holes is the union of region II and the instability region. This stems 
from the fact that the locally stable black hole, at fixed (T, J, Q), is the one with maximum 
entropy, and that fl{S, J, Q) is a decreasing function of 5* at fixed J and Q. 



Unique Black Hole Phase q = o_o5 




FIG. 5. Phase diagrams with Q = 1/1 curve, for charges (5 = and Q = 0.05. The second 
diagram is characteristic for all black holes with < Q < Qc- The dotted curve represents the line 
O = l/l. 

The full phase diagram, taking into account the restriction on the angular velocity, is re- 
ported in figure ^ It implies that part of the small/large black hole phase transition has no 
dual CFT analogue, as in this region one has VL > l/l, and so the Einstein universe in which 
the dual CFT lives would have to rotate faster than light. However the second order critical 
point lies in the Q < l/l region. Furthermore, for charged black holes, the first order phase 
transition occurs aX Q <l/l for sufficiently small angular momentum, and has therefore an 
analogue on the CFT side. 

E. Grand- Canonical Ensemble 

As described in Section |1V B| , the thermodynamics of the grand-canonical ensemble can 
be extracted from the Gibbs potential G{T,Q, $). The natural variables for this ensemble 
are the intensive parameters T, Q and however, to simplify the analysis, we shall study 
the thermodynamics in the (r+, a)-plane, at fixed electric potential $. 
In this plane, the black holes with genuine event horizon are restricted to the region 
$) < a < a+tr('^+5 ^) with a < /, where 

aLrir-,, $) = ^ (r^ - l' + 2/^$^ ± y'^K^Tr + M^) • (63) 

Between these two curves the black holes have a bifurcate horizon, on the curves they are 
extremal, and outside the solution shows a naked singularity. If |<I>| < 1, a^^j, is negative, 
and the entire region below the curve a^tj,(r+, $) represents black hole solutions. For $ > 1, 
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(^extr is positive for < r+ < /a/ — if lies in this range, a~xtr determines a new 
region with naked singularities, and the black holes are found between the two curves. For 
r_|_ outside of this interval, all the points represent black holes, as long as a < Z. Finally, 
in the limit $ — > cxd, the two curves merge and only a line of extremal black holes with 
r+ = a < / survives. 

Another relevant information in the thermodynamic diagram is given by the region where 
the Killing vector dt + ^Hd,j> is timelike all the way out to infinity. It is easy to see that this 
condition is satisfied for a < r^/Z, and so is always fulfilled for r_|_ > /, as we assumed a <l. 
One readily verifies that the condition a < r^// is equivalent to ^7 < 1//. In the diagrams, 
the region at the right of the curve a = r'^/l is the region where the black hole can be in 
thermal and dynamical equilibrium with radiation, and the rotating Einstein universe at the 
boundary spins slower than light. 

Let us now turn to the stability analysis of the black hole. We shall begin with local stability. 
The curve separating stable from unstable regions is given by Eq. (pOf). In fact this equation 
can be solved for a(r+, $) after having substituted q in terms of $ by means of (pOD; the result 
is quite comphcated, but nevertheless adequate for visualizing the local stability regions in 
the {tj^, $)-plane. 

However locally stable black holes can decay into other states with lower free energy. In 
the case of the grand-canonical ensemble, another possible state is given by AdS space, 
which can have arbitrary electric potential $, constant all over the spacetime. The Gibbs 
free energy of AdS space filled with thermal radiation vanishes identically, hence the global 
stability can be investigated by studying the sign of G{T, Q, $). The free energy of the black 
hole is negative for 




HP 



(a,$), (64) 



which means that the black hole dominates in this region. For r+ = rHp(a, $) there is a 
first order phase transition line, with a discontinuity in the entropy, and for < rHp(a, $) 
the AdS solution is globally preferred. However we stress that this does not exhaust the 
question of global stability for these black holes, as other solutions with lower free energy 
may exist ||Tl]]. Indeed, there are regions in the (r_|_, a)-plane in which the black hole, though 



being unstable, has lower free energy that AdS space; it is plausible that a new solution 
exists which minimizes the Gibbs potential. 

Finally we observe that the Q = l/l curve, the local stability curve and the Hawking- Page 
transition curve intersect at the point {r^ = l,a = I) for any choice of the electric potential. 
For $ = $c = 1 these three curves coincide; for $ < $c we find the Q = l/l curve, the 
local stability curve and the Hawking-Page transition curve in this order from the left to the 
right, while for $ > $c they appear in the reverse order. 

We have summarized this analysis in figure |^, which shows the phase structure in the (r_|_, a)- 
plane for various values of the electric potential. We see that for $ < $c = 1? the Hawking- 
Page transition curve is the rightmost one; it separates a large r+ phase of large black 
holes and an AdS phase, which dominates over black holes with small r+. Inside the AdS 
phase region, with increasing r+, we first have a zone of naked singularities separated by 
the extremality line from a region with genuine black holes, which however do not possess a 
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boundary CFT analogue, as they have Vt > l/l, and which are furthermore locally unstable. 
In the next region we find locally unstable black holes which rotate slower than light, and 
finally a region with locally stable but globally unstable black holes which ends on the 
Hawking-Page transition curve. These regions are anyway dominated by the AdS phase. 
For $ > $c the situation is more subtle; a new branch of extremal black holes appears 
and bounds a region of naked singularities for small a. For small r+ we again have the 
AdS dominated phase, but at the right of the Hawking-Page transition curve we find (a) a 
region with locally unstable black holes with f2 > having however negative Gibbs free 
energy, so they are preferred with respect to AdS space, (b) a region with stable black holes 
with VL > l/l, and finally (c) a region of stable black holes with Q < l/l. This raises the 
question for region (a) whether there exists at all a thermodynamically stable state for this 
range of parameters. Furthermore, the region defined by a < a^^-j.(r_|_, $) represents naked 
singularities which have a lower Gibbs free energy than AdS space itself. 



V. CONCLUDING REMARKS 

In the present paper, we studied the thermodynamics of the four-dimensional Kerr- 
Newman-AdS black hole both in the canonical and the grand-canonical ensemble. Thereby 
we encountered an interesting phase structure. In the canonical ensemble, we found a small- 
large black hole first order phase transition, which disappears for sufficiently large electric 
charge or angular momentum. In the points where the small-large transition disappears, 
the first order phase transition degenerates to a higher order one. As pure AdS space is 
not a solution of the Einstein-Maxwell field equations with fixed electric charge, we cannot 
compare the black hole Helmholtz free energy with that of pure AdS space, in order to 
study global stability. It would therefore be desirable to have another reference background 
to our disposal, which describes AdS space filled with charged rotating particles. Another 
possibility would be to have a KNAdS black hole surrounded by a gas of particles carrying 
part of the total charge and angular momentum (cf. the discussion in [|Tl|] ). 

As AdS space can have arbitrary constant electric potential, we can compare the Gibbs 
free energy of the black hole with that of AdS space in the grand-canonical ensemble. In 
this way, we found the Hawking-Page transition curve. Also in this context, some intriguing 
features arise, e. g. one finds black holes which are locally unstable, but nevertheless have 
lower Gibbs free energy than AdS space itself. There is also a region of naked singularities 
in the (r_|_, a)-phase diagram, which are energetically preferred with respect to AdS space. 
It would be interesting to study if these naked singularities play some role in the dual CFT. 

We hope to have clarified that a consistent treatment of rotating AdS black hole ther- 
modynamics requires the usage of the angular velocity Q of the rotating Einstein universe, 
in which the dual conformal field theory lives. This angular velocity differs from that of the 
event horizon, which enters the thermodynamical description of asymptotically flat black 
holes. The reason for this discrepancy is that the angular velocity (|15D appearing in the 



canonical form of the metric, does not vanish at inflnity. (It even has the opposite sign of 
Qh, cf. (|15D ). The fact that a consistent thermodynamics can be formulated using fl is in 
full agreement with the AdS/CFT correspondence. 

The black holes we have considered can also be lifted to eleven dimensions, as solutions of 
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= 1, D = 11 supergravity |[T0|j39[] . These solutions represent the near- horizon hmit of M2- 
branes with spherical worldvolumes, rotating both about an axis of the worldvolume and in 
planes orthogonal to the worldvolume, the latter rotation giving rise to four R-charges upon 
Kaluza-Klein compactification. The Kerr-Newman black holes considered here correspond 
to the special case where all transverse angular momenta of the M2-brane, and hence all 
R-charges of the black hole, are equal 0. 

The thermodynamic structure we have elaborated, corresponds then to that of the world- 
volume field theory of coincident such branes. 

It would be very interesting to give a microscopic description of the Kerr-Newman- AdS 
black holes in terms of this dual CFT. Let us recall in this context that only black holes 
whose angular velocity satisfies 

Q<j (65) 

correspond to thermal states in the dual field theory, as otherwise the rotating Einstein 
universe, in which the dual CFT lives, would rotate faster than light. One easily shows that 
for zero temperature black holes, condition (]65|) implies 



> a{a + ly (66) 



for the electric charge and rotation parameter appearing in the metric. If the bound 
is saturated, we have a BPS state. It would be interesting to study if the corresponding 
field theory at weak coupling reproduces the Bekenstein-Hawking entropy of the black hole 
in this BPS limit. As the M2-brane worldvolume field theory is the strong coupling limit 
(conformal point) of the D2-brane worldvolume theory the relevant gauge theory would 
be that on the worldvolume of D2-branes. In any case, in the presence of a global background 
L'^(l) current (corresponding to the black hole charge Q), the zero temperature state of the 
dual field theory must be highly degenerate, as the extremal black hole has nonvanishing 
entropy^. It would also be interesting to consider the five-dimensional case, corresponding 
to A/" = 4, D = 4 SYM on a rotating Einstein universe |T3]. However, the generalization of 



the rotating five- dimensional AdS black holes found in |jT8| to nonzero charge seems to be 
a nontrivial task. If such a solution were known, one could study if there exist BPS states 
(which can only occur for nonzero charge |l^), and compare their entropy with the one 
obtained on the CFT side. Apart from the possiblity of identifying black hole microstates, 
this would also be a check of the AdS/CFT correspondence. 



thermodynamical discussion of branes rotating in planes orthogonal to the worldvolume can 
be found in ||40|- 43 , 191] . In [40|, a field theory model of rotating D3-branes was proposed. This was 
later generalized in to rotating M-branes. 



^Note that for zero U{1) current, the extremal black holes have Q > l/l [jT8|, and therefore do 
not correspond to states in the CFT. 
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Grand-canonical ensemble, <!> = 



Grand-canonical ensemble, (1> = 0.5 




FIG. 6. Phase structure in the grand-canonical ensemble in the (r+, o)-plane, for various values of 
the electric potential For $ = we have, from the left to the right, the Q = l/l curve, the Une 
bounding the stabihty region, and the Hawking-Page transition curve. For < $ < $c = 1, in 
addition to these three hnes, the extremahty curve appears at the left. For the critical potential 
$ = <I>C) the Q = l/l curve, the stability bound, and the Hawking-Page transition line coalesce. 
Finally, for <1> > <I>c, an additional extremality line appears (the one intersecting the r+ axis twice), 
and the order of the Q = l/l curve, the stability bound, and the Hawking-Page transition line is 
reversed. The various regions bounded by these curves are discussed in the text. 
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